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Abstract. Observational evidence of a hydrodynamically evaporating upper atmosphere of HD209458b (Vidal-Madjar et al. 
2003; 2004) and recent theoretical studies on evaporation scenarios of "Hot Jupiters" in orbits around solar-like stars with the 
age of the Sun indicate that the upper atmospheres of short-periodic exoplanets experience hydrodynamic blow-off conditions 
resulting in loss rates of the order of about 10 10 -10 12 g s (Lammer et al. 2003; Yelle 2004; Baraffe et al. 2004; Lecavlier 
des Etangs et al. 2004; Jaritz et al. 2005, Tian et al. 2005; Penz et al. 2007). By studying the effect of the Roche lobe on the 
atmospheric loss from short-periodic gas giants we found, that the effect of the Roche lobe can enhance the hydrodynamic 
evaporation from HD209458b by about 2 and from OGLE-TR-56b by about 2.5 times. For similar exoplanets which are closer 
to their host star than OGLE-TR-56b, the enhancement of the mass loss can be even larger. Moreover, we show that the effect 
of the Roche lobe raises the possibility that "Hot Jupiters" can reach blow-off conditions at temperatures which are less than 
expected (< 10000 K) due to the stellar X-ray and EUV (XUV) heating. 
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IVidal Madi ar" st al. I d2003). different blow-off scenarios for evaporating hydrogen-rich atmospheres are discusse d in the literature 
dSasselov 120031: iLammer et al. feOollLecavelier des Etangs et al. l2004tlYelle l2004HVidal Madiar et al. l2004tlGrieBmeier et al. 



(Sasselov ZUU3; Lammer et al. ZUUJ; Lecavelier des Utangs et al. ZUU4; Yelle /UU 
20041 iBaraffe et al. 112004 lErkaev et al. ll2005blTian et al. ll2005HPenz et al. 1 12006). 



Because hydrogen-rich upper atmospheres of short-periodic gas giants in orbits with semi-major axes less than < 0.1 AU can 
be heated to temperatures of up to 10000-20000 K by X-rays and EU V (XUV) radiation of the central star (Lammer et al. 2003; 
Yelle 2004), they experience hydrodynamic blow-off. lYelle (2004) studied in detail the photochemistry of the thermospheres 



of "Hot Jupiters" and found that the lower thermosphere may be cooled primarily by infra r ed radiative emissions from Ht 
molecular ions which are cre ated by photoionization of H2 molecules and related ion chemistry (Yelle l Eool) . 



Rec ently, iLecavelier des Etangs et al." I d2004 and Jaritz et al. (2005) discussed the Roche lobe effects on hydrodynamic mass 



loss from the "Hot Jupiters" atmospheres and argued that for some close-in exoplanets, due to the expected high exospheric 
temperatures, the exobase level r exo can reach the Roche lobe before classical hydrodynamic blow-off conditions may develop 
and, thus, affect their atmospheric loss rates. The Roche lobe is defined as the last equipotential surface around a planet where its 
gravitational potential energy is zero. Below this boundary, the surfaces of c onstant potential encompass the ho st star. Because 



atmospheric particles at and beyond the Roche lobe can escape unhampered. ILecavelier des Etangs et al. (2004) argued that the 



Roche lobe can be seen as an equivalent to the exobase level. However, there is an important difference between the Roche lobe 
and the exobase, since at the Roche lobe and above it a planetary atmosphere is not bound to its planet by gravitational forces 
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and can freely escape, whereas at the exobase level (where r exo < t-ri) gas temperature should be above some critical value for 
blow-off to occur. 

In the typical blow-off scenarios considered, for example, by Watson et al. (1981) and Jaritz et al. (2005) it is assumed that the 
stellar X-ray and EUV radiation is absorbed mainly at the so-called expansion radius r\ (Watson et al., 1981) or at r exp (Jaritz et 
al., 2005) which is located for most "Hot Jupiters" at a planetocentric distance close to r R ] or above. Therefore, Jaritz et al. (2005) 
proposed that, depending on a planetary and stellar mass, planetary size, stellar type and orbital distance, both scenarios, the 
classical hydrodynamic blow-off, and the so-called "geometrical blow-off", as expected by Lecavelier des Etangs et al. (2004), 
may occur. 

However, hydrodynamic model simulations by Yelle (2004) which include photochemistry of a hydrogen-rich atmosphere, 
indicate that the optical density for the stellar XUV radiation absorption Txuv by evaporating hydrogen is <k 1 at t-ri or at 
distances corresponding to Watson's r\, so that the main part of the XUV radiation is absorbed at lower altitudes. This is also 
in agreement with a recent study by Kulikov et al. (2006) who showed that in a terrestrial type planetary atmosphere even at 
very high XUV fluxes the bulk of the XUV radiation is absorbed at a distance rxuv (where Txuv = 1) which is much closer 
to the planetary radius r p \ than r\. At altitudes which are much higher than r p i or r X uv> the absorption optical depth t X uv of 
an evaporating atmosphere is much less than 1. And below r X \jv, where the optical depth Txuv » 1, in situ XUV atmospheric 
heating becomes negligible. 

The aim of this paper is to show that for some of the observed close-in "Hot Jupiters' the Roche lobe can be fairly close to a 
planet and its exobase level, which can substantially enhance the mass loss of a planetary atmosphere as compared to the classic 
blow-off conditions, and may thus have important evolutionary implications. 



2. Roche lobes of close-in exoplanets 



We consider two spherical masses, M p \ (exoplanet) and M star (star), separated by an orbital distance d, which rotate around their 
common center of mass. In the rotating coordinate system, the energy per unit mass of a test particle in the ecliptic plane is given 
by Paczyhski (1971) 



GM pl GM star G(M pl +M star )s 2 



2d 3 



(1) 



where s is the distance from a given point to the center of mass, r a and r\, are distances to the centers of the exoplanet and the 
star, respectively, and G is Newton's gravitational constant. 

The first term in eq. (1) represents the potential of the exoplanet and the second the potential of the star. The third term is the 
result of the orbital motion of the whole system. By introducing dimensionless quantities 5 — M p i/M star , A - d/r p \, r\ - r a /r p \, 
we analyze the variation of the potential along the axis which connects the exoplanet with its host star 
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For this potential, there are two locations, the Lagrangian po int LI and the L agrangian point L2 which are saddle points. Both 
points are rather close to each other for a small ratio of mass ( Gu et al.ll2003 ) 



7"L1,L2 



d, 



therefore, we consider as the Roche lobe boundary t-ri, approximately 
r Rl * (|) ' d. 

The potential difference between the Roche lobe boundary and the planetary radius (surface) can be written as 



(4) 



(5) 



AO = Oo 
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Exoplanet 


star type 


Mi, [M ] 


M p i [Mj„ p ] 


''pi [Hup] 


rm [/pi] 


l/K 


K 


d [AU] 


HD209458b 


GOV 


1.05 


0.69 


1.43 


3.9 


2.10 


0.47 


0.045 


OGLE-TR-56b 


G 


1.04 


1.45 


1.23 


3.0 


2.39 


0.41 


0.023 


OGLE-TR-132b 


F 


1.34 


1.01 


1.15 


3.5 


1.90 


0.52 


0.031 


OGLE-TR-113b 


K 


0.77 


1.35 


1.08 


3.6 


1.75 


0.57 


0.023 


TreS-1 


KOV 


0.87 


0.75 


1.08 


4.9 


1.48 


0.67 


0.023 


OGLE-TR-lllb 


GorK 


0.82 


0.53 


1.00 


5.8 


1.34 


0.74 


0.047 


OGLE-TR-lOb 


GorK 


1.22 


0.57 


1.24 


3.7 


1.9 


0.52 


0.042 



Table 1. Factors l/K and K for 7 "Hot Jupiters" where t he planetary mass and radius is known. T he planeta ry and s tellar 
param eters for HD209458 b dBarman et al l2002t lVid al Madiar et al. 120031). OGLE-TR-5 6 b dBurrows et al. l2004lBaraffe et al I 
120041) . OGLE-TR -132 b dMoutouetal Il2004. OGLE-TR-113 b dBouchv et al. Il2004. OGL E-TR-111 b dPontetal.ll2004. 
OGLE -TR-132 b dBouchv et al. Il2004 iMoutou et al. Il2004 . OGLE-TR-10 b dBouchv et al. 112005b and TreS-1 dAlonsoetal " 
2004) were used for the calculation of rm and K. 



By assuming d » tri > r v \ and M star » M v \ what means X » tj > 1 and 6 <k 1, we can simplify the expression for the potential 
difference 



AO = O 



x _ 3 77(1 + 77) 



2 (M 3 



Introducing the dimensionless quantities A and 77 in eq. (5), we can further simplify eq. (7) for the potential difference to 
(77 - 1) 2 (2t7 + 1) 



AO = O, 



0- 



2?7 3 



(7) 
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3. Implications for atmospheric blow-off 

By considering a hydrodynamic regime of atmospheric escape we apply the energy conservation equation 



mv 5 
mAO+ — + -k(T m 



T ) 



f 



qr 2 dr, 



(9) 



where F is the loss rate of particles per steradian, q = gxuv ~<?ir, with the XUV volume heating rate gxuv and the cooling rate 
due to IR emitting molecules like (Yelle, 2004), m is the mass of the evaporating particles, and v is the outflow bulk velocity 
at the Roche lobe boundary r^, k is the Boltzmann constant, Tri is the temperature of hydrogen at the Roche lobe and Tq is 
the temperature at a distance close to r p i, which is about the effective radiative temperature of the exoplanet. Substituting the 
potential difference AO into (9), we obtain 



r = 



rqr 2 dr 



'pi 



+ 2f + U(T m -T ) 



where the factor K is 

R= Iftfr + 1) 
(2/7 3 ) 



(10) 



(11) 



Neglecting the kinetic energy term mv 12 and also the thermal energy 5k/2(Tm - To) and introducing the planetocentric distance 

n 



( 



''1 



r qrdr 



(12) 



we can obtain the energy limited escape rate equation similar to that derived by Watson et al. (1981), which corresponds to the 
distance below which the incoming stellar XUV radiation 7xuv is absorbed by an evaporating atmosphere 



r = 



fpl^/xuv 

mM p \GK' 



(13) 



Fig. 1 illustrates the difference between the XUV absorbtion radius r\ of Watson et al. (1981), the radius rxuv where the bulk 
XUV flux is absorbed and the energy absorption function qxuv has its peak, and the Roche lobe r&\. Depending on the planetary 



4 



N. V. Erkaev et al.: Roche lobe effects on the atmospheric loss of "Hot Jupiters" 




Fig. 1. Illustration which compares r p \, rxuv, f^a, <7xuv with the XUV absorbtion radius r\. If the exoplanet experiences hydrody- 
namic conditions, blow-off occurs at the altitude level rxuv where the XUV energy absorption function qxuv has its maximum. 
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0.1 I ■ ■ ■ < i i i I 
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Fig. 2. Roche lobe induced factor K and atmospheric mass loss enhancement l/K as a function of rm normalized to r v \. 



and stellar parameters, r\ can be located outside the Roche lobe. Due to the low number density of the evaporating hydrogen 
at planetary distances around and above tri, the optical density of the XUV radiation txuv is <k 1 and the main part of the 
stellar XUV radiation is absorbed at r X uv which is close to r v \. This altitude r X uv corresponds generally to the thermosphere and 
ionosphere. 

If the bulk of the stellar XUV radiation is absorbed in the thermosphere at altitudes r X uv < r m, the particle loss can be 
enhanced substantially by the effect of the Roche lobe which manifests itself in the loss enhancement factor 1/K in eq. (13). 
Table 1 and Fig. 2 show the Roche lobe induced mass loss enhancement factor 1 /K (and also K) as a function of tri in planetary 
radii r p \. 

One can see from Table 1 that OGLE-TR-56b at an orbital distance of about 0.023 AU experiences the strongest enhancement 
of the Roche lobe affected evaporation resulting in a factor of about 2.4, while the Roche lobe induced enhancement of mass loss 
at HD209458b at 0.045 AU is about 2. As one can see from Fig. 2, the mass loss of "Hot Jupiters" which orbit at < 0.02 AU 
around their host stars, could be dramatically larger, because tri would move closer to r p \ and that enhances the evaporation. One 
more important effect due to the Roche lobe on "Hot Jupiters", which experience hydrodynamic conditions, manifests itself in 
the potential energy difference which is less than the thermal energy at the exobase for blow-off to occur 

mAO < kT m , x , (14) 
Here r max is the maximum exobase temperature produced by the XUV heating. This equation yields 
GmM v \K 



< I- (15) 



One can see that the effect of the Roche lobe helps to satisfy this condition for lower temperatures than expected for the classic 
blow-off (that is for K — 1). In other words, if the blow-off temperature for an exoplanet without the effect of the Roche lobe 
(K x I) is, for example, 10000 K, a similar exoplanet, but which is closer to its host star, may start to evaporate hydrodynamically 
due to the Roche lobe effect at about 5000 K if the factor K is « 0.5 like it is for OGLE-TR-132b (see Table 1). 

This result is very important, because it indicates that the effect of the Roche lobe can enhance the possibility that "Hot 
Jupiters" may reach hydrodynamic blow-off conditions more easily, even if their atmospheres have a high amount of molecules 
like Ho , which act as IR-coolers in the thermosphere. Both effects, the enhanced mass loss and the higher probability that "Hot 
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Jupiters" reach hydrodynamic blow-off conditions at very close orbital distances to their host stars, may enhance the evaporation 
rate. The results of our study have to be included in the statistical mass-radius analysis of hot exoplanets expected to be detected 
during the CoRoT mission in the near future. 

4. Conclusions 

Our study shows that hydrodynamically driven atmospheric mass loss from "Hot Jupiters" at close orbital distances that are 
much less than 0.05 AU, may be strongly enhanced due to the Roche lobe effect as compared to non-affected exoplanets. We 
found that the mass loss due to the Roche lobe can be enhanced several times if the Roche lobe is located closer to a planet 
at a distance of a few planetary radii. Furthermore, our study indicates that the Roche lobe effect may also help "Hot Jupiters" 
to attain hydrodynamic blow-off conditions even if their exospheric temperatures are lower than those required for the blow-off 
to occur in the case of a classic Newtonian gravitational potential of a planet. Both effects may have a strong impact on the 
atmospheric evolution of short periodic hydrogen-rich gas giants. 
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